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Abstract 

Using the basic ingredient of supersymmetry, we develop a simple alter- 
native approach to perturbation theory in one-dimensional non-relativistic 
quantum mechanics. The formulae for the energy shifts and wave functions 
do not involve tedious calculations which appear in the available perturba- 
tion theories. The model applicable in the same form to both the ground 
state and excited bound states, unlike the recently introduced supersym- 
metric perturbation technique which, together with other approaches based 
on logarithmic perturbation theory, are involved within the more general 
framework of the present formalism. 

1 Introduction 

An exact solution of the Schrodinger equation exists only for a few idealized 
problems; normally it has to be solved using an approximation method such as 
the perturbation theory (PT), which constitutes one of the most powerful tools 
available in the study of quantum mechanics in the atoms and molecules. PT 
is applied to those cases in which the real system can be described by a small 
change in an exactly solvable idealized system. In this form we can describe 
a great number of problems encountered especially in atomic physics, in which 
the nucleus provides the strong central potential for the electrons; further inter- 
actions of less strength are described by the perturbation. Examples of these 
additional interactions are: the magnetic interaction (spin-orbit coupling), the 
electrostatic repulsion of electrons and the influence of external fields. But in 
spite of widespread application of this theory, its basic analytical properties are 
poorly understood. One of our objectives in this paper is to illustrate selected 
important aspects of the perturbation theory within the frame of supersymmetric 
quantum mechanics. 

Performing explicit calculations in non-relativistic quantum mechanics using 
the familiar Rayleigh-Schrodinger perturbation expansion is rendered difficult by 



the presence of summations over all intermediate unperturbed eigenstates. Al- 
ternative perturbation procedures have been proposed to avoid this difficulty, 
notably the logaritmic perturbation theory (LPT) [l]-[4] and the Dalgarno-Lewis 
technique [5]- [8]. The virtue of LPT is its avoidance of the cumbersome summa- 
tion over states for second- and higher-order corrections in Rayleigh-Schrodinger 
perturbation theory. Unfortunately, it has problems of its own in calculating 
corrections to excited states, owing to presence of nodes in the wave functions. 
Various schemes have been proposed to circumvent the resulting singularities 
[4, 9, 10]. 

Such is the status of LPT after over 20 years of active development. Mean- 
while, supersymmetric quantum mechanics (SSQM) [11, 12] has developed im- 
mensely since the first models were introduced [13, 14]. Several approximation 
methods using SSQM formalism have been developed, including the supersym- 
metric perturbation theory (SSPT) of Cooper and Roy [15]. Recently, Lee [16] 
has shown that SSPT and LPT are entirely equivalent and fortuitously, each 
turns out to resolve difficulties encountered in the other. Namely, LPT formulas 
for energy corrections obviate tedious procedures in the SSQM method, while 
the use of SSQM partner potentials with virtually identical bound state spectra 
solves difficulties with excited states encountered in LPT. Although the iterative 
procedure in SSPT may not actually reduce the calculational workload, it does 
cast the calculations into a physically-motivated, visualizable framework. 

In this letter, starting from the first priciples, we develop a more economi- 
cal scheme which yields simple but closed perturbation theory formulae leading 
to the Riccati equation from which one can actually obtain all the perturbation 
corrections to both energy level shifts and wave functions for all states, unlike 
the other models mentioned above. In the application of the present method to 
the n th excited state, one requires knowledge of the unperturbed eigenfunction 
Xn(f) but no knowledge of the other eigenvalues or eigenfunctions is necessary. 
The procedure introduced here does not involve either tedious explicit factoring 
out of the zeros of x«( r ) [1> 2] or introduction of ghost states [4] as were the 
cases encountered for applying LPT to excited states. Since, the present model 
offers the explicit expressions for the energy corrections, which are absent in the 
original SSPT while the treatment of Lee [16] for such calculations has mathe- 
matical complexity, and provides a clean route to the excited states, which are 
cumbersome to analyze in both LPT and SSPT, our results can be thought of as 
a generalization of logarithmic and supersymmetric based perturbation theories. 
This is the another objective in the present work. 

In the following section we introduce the model and discuss briefly the phyiscs 
behind the formulation. In Section 3, some applications are given and the power 
of the present technique is illustrated when compared to the calculation technique 
of other theories considering the whole of states. Some concluding remarks and 
summary of the work are drawn in the last section. 
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2 Formulation 



The goal in SSQM is to solve the Riccati equation, 

h 



W 2 (r) 



'2m 



W\r) = V(r) - E , 



(1) 



where V(r) is the potential of interest and E Q is the corresponding ground state 
energy. If we find W(r), the so called superpotential, we have of course found 
the ground state wave function via, 



ipo(r) = N exp 



2m 



h 



W(z)dz 



(2) 



where N is the normalization constant. If V(r) is a shape invariant potential, we 
can in fact obtain the entire spectrum of bound state energies and wave functions 
via ladder operators [11]. Through the article, this basic ingredient of SSQM given 
by (1) and (2) will be extended and used for the treatment of excited states. 

Now, suppose that we are interested in a potential for which we do not know 
W(r) exactly, and the corresponding Hamiltonian is not factorizable but almost 
factorizable. More specifically, we assume that V(r) differs by a small amount 
from a potential V (r) plus angular momentum barrier if any, for which one solves 
the Riccati equation explicitly. For the consideration of spherically symmetric 
potentials, the corresponding Schrodinger equation for the radial wave function 
has the form 



k 2 m 

2m ifj n {r) 



[V{r) - E n 



V(r) 



V (r) + 



h 2 



2m 



+ AV(r), (3) 



where AV is a perturbing potential. Let us write the wave function ip n as 

*l>n(r) = Xn(r)(j)n(r) , 



(4) 



in which Xn is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas n is a moderating function corresponding to the perturbing 
potential. Substituting (4) into (3) yields 



2m \Xn <P n Xn rn / 



(5) 



Instead of setting the functions Xn and <p n , we will set their logarithmic derivatives 
using the spirit of Eqs. (1) and (2); 



h Xn 

/2m Xn 



/ 2m(/)n 



(6) 
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which leads to 



T) 1 v" h 



2m Xn ' V2m 



V (r) + 
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2m r 



2 



(7) 



where e n is the eigenvalue of the unperturbed and exactly solvable potential, and 
+ = AW 2 n - ~^=AW' n + 2W n AW n = AV(r) - Ae n , (8) 

2m \(p n Xn (pnj \j2m 

in which Ae n is the eigenvalue for the perturbed potential, and E n = e n + Ae n . 
Then, Eq. (5), and subsequently Eq. (3), reduces to 

{W n + AW n ) 2 - -A= (W n + AW n )' = V - E n , (9) 
Vim 



which is similar to Eq. (1). In priciple as one knows explicitly the solution of 
Eq. (7), namely the whole spectrum and corresponding eigenfunctions of the 
unperturbed interaction potential, the goal here is to solve only Eq. (8), which 
is the main result of this letter, leading to the solution of Eqs. (3) and (9). 

Eq. (8) is a closed analytical form in comparing to lengthy SSPT and LPT 
expressions, in particular for the excited states. In this respect, the present 
formulation has a more general form than the available perturbation theories. 
Though this point will be clear in the next section through the applications, it 
would be convenient at this stage to clarify how Eq. (8) involves in a compact 
form the supersymmetric and logarithmic perturbation theory expressions. As 
the equivalance of SSPT to LPT has already been clarifed [16], we consider here 
only the framework of SSPT and show that SSPT is a subset of the present model. 

For the perturbation technique, we have initially assumed that we could split 
the given potential in two parts, Eq. (3). The main part corresponds to a shape 
invariant potential, Eq. (7), for which the superpotential is known analytically 
and the remaining part is treated as a perturbation, Eq. (8). If necessary, one 
can expand the functions related to the perturbation in terms of the perturbation 
parameter A, 

oo 

AV(r;X) = ]TA fc A\4(r), 
fc=i 

oo 

AW n {r-\) = Y, xkAW nk(r) , 



fc=i 

oo 



Ae n (A) = ]TA fc e„ fc , (10) 

k=l 

where A will eventually be set equal to one. Substitution of the above expansion 
into Eq. (8) by equating terms with the same power of A on both sides yields up 
to O (A 3 ) 

2W n AW nl - ~^=AW' nl = AV 1 - Ae nl , (11) 



'2m 
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AW 2 , + 2W n AW n2 - -=A< 2 = AV 2 - Ae n2 , (12) 

V2m 

2 (W n A^ n3 + AW nl AW n2 ) - -^AW; = A\/ 3 - Ae n3 , (13) 

V2m 

which are exactly SSPT expressions appeared in [15]- [18] but for the case n — 0. 
Eq. (8) and its expansion, Eqs. (11-13), give a flexibility for the easy calcula- 
tions of the perturbative corrections to energy and wave functions for the n th 
state of interest through an appropriately chosen perturbed superpotential, un- 
like the other perturbation theories. We will show in the next section through the 
applications that this feature of the present model leads to a simple framework 
in obtaining the corrections to all states without using complicated and tedious 
mathematical procedures. 



3 Application 

3.1 perturbed harmonic oscillator 

Let us start with the most elementary kind of perturbation calculation. Consider 
a perturbed harmonic oscillator potential of the form 



:i ill -- ■ Li; 1 " 

V(r;X) 



1 22 £(£+ l)h 2 
■mw r + 



+ -m\w 2 r 2 , (14) 



2 2mr 2 

which woud arise by increasing the spring constant of a harmonic oscillator from 



K to (1 + \)K since w = yK/m. The perturbed potential in its present form 
can of course be solved exactly. Neverthless, we want to solve this problem in 
the light of our approach to test the effectiveness of the present technique. This 
application will also clarify that our derivation is much simpler and more direct 
than the other methods, and hence provides a useful alternative. 

The whole spectrum and corresponding wave functions for the unperturbed 
part of the problem is well known in the literature. Then, starting with the nor- 
malized wave function of the harmonic oscillator system for any state of interest, 
one can easily calculate the corrections at one step to the energy and eigenvalues 
for the perturbed potential setting an appropriate superpotential satisfying Eq. 
(8). For the comparison of our results with those of recent works, we will perform 
the calculations here for the ground state n — 0. Starting with Eq. (8), 

Ti 1 
AW 2 =0 - — =A< =0 + 2W n=0 AW n=0 = -m\w 2 r 2 - Ae n=0 , (15) 
V2m I 

one readily sees that the perturbed superpotential (AW) satisfying the above 
equation has the form 

AW n=0 = ^wr (VTTX - l) , (16) 
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since, from the literature [11], the superpotential W corresponding to the unper- 
turbed potential is in the form 

(17) 

Upon substituting (16) and (17) into (15) yields 

Ae n=0 = (t+l) (v/TTA - l) hw, (18) 
and using (16) with (2) leads to the moderating function for the ground state 



n =o = exp 



mw 2 



2h 



r 2 (l - VT+X) 



(19) 



As the ground state wave function Xn=o corresponding to the superpotential in 
(17) and energy e n=0 for the normalized unperturbed harmonic oscillator are 



at /+i ( rnwr 2 \ ( 3\ , 

Xn=o = N r + exp I — I , e n=0 = 1 1 + - J hw , 



(20) 



the total wave function ■0n=o f° r the whole system in (14) and corresponding 
energy are obtained as 



V>n=0 = Xn=O0n=O = N / +1 eX P ^ + ^ , E n =0 = V 1 + A (V + ^ ^ . 

(21) 

These are indeed the correct results which can be verified explicitly by the use of 
Eq. (9), since this problem has an exact analytical solution. 

Lee [16, 17] recently has studied this problem within the frame of supersym- 
metric perturbation theory and worked out it up to the second order. From the 
expansion of the closed expressions in (21) in A such that y/1 + A = 1 + A/2 — 
A 2 /8 + . . ., it can clearly be seen that Lee's treatment indeed appears as a subset 
of the present calculation procedure. In addition, our sophisticated approach nei- 
ther involves cumbersome procedures nor tedious calculations. In what follows, 
we further show that the present model is also applicable in the same form to 
bound excited states without any difficulty, unlike the works in [15]-[17]. 

As the potential considered in (14) is a shape invariant potential, using (9) 
together with (16) and (17) for n = and having in mind the shape invariance 
property for the exactly solvable potentials [11], 

V+(r, oq) = V-(r, cii) + i2(ai) ; R (a x ) = 2/Wl + A , (22) 

where V_ is the exactly solvable potential in (9) and V + is its supersymmetric 
partner 



(W n=0 + AW n=0 ) 2 + (W n=0 + AW n=0 )' = V+- E n=0 , (23) 

V2m 
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and a = £, ai = t + 1. Therefore, 

£ n = 2 i?(a s ) + (i + hw = (2n + £ + ^ kwVTT\ , (24) 

which are the whole spectrum for the perturbed harmonic potential in (14). This 
very simple procedure removes difficulties encountered in SSPT and LPT in deal- 
ing with excited states. Again, this result agrees with Eq. (47) of [16] in which 
the excited state energies are presented up to only second order due to nasty 
calculation procedure. 

To find the energy of the excited states, all we had to do was to perform 
calculations using a proper superpotential corresponding to the ground state wave 
function of the perturbed potential. The LPT literature goes to great pains to find 
ways to avoid using excited state wave functions. The study of present example, 
in particular the procedure used through Eqs. (22) to (24), thus illustrates in 
this respect how SSQM allows us to use the LPT based formulas for the entire 
spectrum, by moving to a partner potential whose ground state coincides with the 
excited state of the potential initillay considered. The reader is referred to [11] for 
a further discussion on supersymmetric partner potentials and their applications. 



3.2 perturbed Coulomb interaction 

Kim and Sukhatme [9] obtained a set of expressions for ground and excited state 
wave functions and energies in perturbation theory (henceforth referred to as 
KS), that do not involve infinite sums, and which they consider a generaliza- 
tion of LPT. Subsequently some expressions in this formalism were simplified 
by Mavromatis [19]. Later, KS approach were carefully compared [20] with the 
results of the LPT formalism and step by step connection between the two for- 
malisms was shown and commented on via an illustrative example, together with 
the discussion on the relation between KS approach to Dalgarno-Lewis formalism 
[5]- [8]. 

As we are aware of the relation between LPT and SSPT, one then may expect 
also the connection between SSPT and LPT based KS approach. Thus, we focus 
through this section and the other examples in the next sections on constructing 
a bridge between the present generalized supersymmetric perturbation formalism 
and KS approach and show explicitly the equivalence between the two aproaches, 
considering of course the expanded form of Eq. (8). 

We proceed first with the present treatment of the simple example used in 
[20], which is the perturbed Coulomb system, 



V(r) 



2mr 2 



e 2 



+ A-, (25) 



where, for a particular example, A denotes the increse in the charge of the nucleus 
in case of an interaction between electron and nucleus in hydrogen like atoms. 
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The superpotential, wave function and energy for the ground state corre- 
sponding to the unperturbed potential, respectively, are [11] 



n=0 



m e 2 (e+l)h 



j+i < me 



2 



Xn=o = Nr exp ['^^r j , 

€n =° ~ ~2h 2 (l + lf (26) 

From Eq. (8), AW n=0 involving all the corrections corresponding to the perturb- 
ing potential AV = Ae 2 /2r is readily obtained as 

AW - = "/I (7TT)i ' < 27 > 

and the the total correction to the enery Ae n= o, together with the moderating 
function due to the perturbing potential are 

. me 4 (\ 2 _\ , / Xme 2 \ 

Ae„. = - Xj , = exp [^^r j ■ (28) 

Hence, the full ground state wave function and energy spectrum for the perturbed 
Coulomb system in (25) are 



^n=0 = Xn=O0n=O = N / +1 exp 



me ( A 



(e + i)h 2 \ 2 



E n - = - ™ A 2 fl--V • (29) 
2ft 2 (£ + l) 2 V 2 / 

One can justify the result above via Eq. (9) and the full excited energy spectrum 
can be easily found using the procedure given through Eqs (22-24), for which the 
superpotential should have the final form W n= o + AW n=0 . 

Now, to see the close relation between the formalism presented in this letter 
and that of Kim and Sukhatme [9], the reader should go back to Eqs. (11-13) 
which are the expansion of the present formalism yielding individual corrections 
in order, having in mind that (p^ s is expanded in orders of A in KS approach. 
In the light of recent works [16, 17] regarding SSPT, one sees after some algebra 
that the expansion terms of moderating function <p^ s in KS technique are related 
to those of the superpotential AW n appearing in (11) through (13). For clarity 
consider only the first order, then 

A ^--i^ s • (30) 
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since, from (11), 

i r 

AW nl (r) = \ / xl(z) [Ae nl - m{z)] dz . (31) 

n Xn{T) -to 

Mavromatis [20], using KS and LPT formalism and working within Ti = m = 
e 2 = 1 unit system for £ = case, investigated the same problem for the ground 
state solution {n = 0) for which after some tedious integrals he arrived at 

K S = , = ^£a 2 , = i! ^A 3 , . . . (32) 

and Ae i = A/2, Ae 2 = — A 2 /8 with higher order energy corrections being 
zero. We should remark at this point that the present approach has provided 
exact result in a closed analytical form, Eq. (28), without dealing with nasty 
integrals, since the moderating function in (32) eventually can be written as 
= exp [A/2 (r — a)] with an extra constant, exp (— Aa/2), though the con- 
stants do not affect the energy, where a comes from the lower limit of the inte- 
grals carried out. Considering Eqs. (30,32) within the same frame and the whole 
discussion given in this section, together with the physics behind Eqs. (8-13), one 
can see that the model introduced here through Eq. (8) unifies not only SSPT 
and LPT but also the KS approach. In the following section, we will make more 
clear this point with the further examples used in [9]. 

3.3 harmonic oscillator with linear perturbation (n=l state) 

Kim and his co-worker considered a one dimensional perturbed harmonic oscilla- 
tor potential with no angular momentum barrier 

V = ^w 2 r 2 + Xr + B , (33) 

involving a linear perturbation term. Picking the particular case of an applied 
uniform electric field to a charged particle moving in a simple harmonic oscillator 
potential provides a simple physical interpretation to the linear term in (33) where 
in this case A stands for the strenght of the field applied. This explanation is also 
valid for the next example considered in the following section but for a different 
unperturbed potential. 

In [9] , the energy level shifts and corrections to the wave function were worked 
out up to the second order for the first excited state {n = 1), 
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For a clear comparison we work within H = 2m = 1 unit system as in [9]. From 
(30) and (34), AW in the first order should read 

AW U = - (l + -V) > ( 35 ) 



w \ wr 2 / 

and the use of (35) in (11) yields explicitly the first order energy shift Ae n = B, 
since 

/ \ IDT" ' 

Xn=i = Nr exp (-wr 2 /4) , W n=1 = - , e n=1 = 3w/2 . (36) 

\ / 2 r 

The second order calculations can be carried out in the same manner for com- 
parison. However, our aim here is to clarify the effectiveness of the technique 
introduced in this article by calculating all corrections explicitly in a simple way. 
For this reason we turn back our attention to the shifted harmonic oscillator po- 
tential in (33) in order to remind the reader that the potential in (33) can be 
exactly solved in its present form by setting the superpotential 

HIT A 

WZ a f = -«- + -> (37) 

from which one sees that B = —w and E^ff 1 = 3w/2 + B — X 2 /w 2 = e n= i + 
Aen + Aei2 which justifies the results obtained above. This makes clear that the 
energy level shifts due to higher orders in A vanish. Eq. (37) leads to the full 
wave function in the light of Eq. (6) 

^exact = Ni exp (_ wr 2/4^ exp (_Ar/u?) . (38) 

Further, it is not difficult to see that Eq. (9) provides the exact full energy 
spectrum involving all excited states using the shape invariance property of the 
potential in (33) without carrying out tedious integrals as in the other pertur- 
bation theories. The physics behind this exact solution would shed a light to be 
able to use Eq. (8) for an easy calculation of the contributions comes from due 
to the perturbed terms, which are briefly discussed below. 

To illustrate once more the power and elegancy of the present formalism we 
here show that all the corrections to the sytem of interest can be calculated 
simply by Eq. (8). For the choice of correct superpotential leading to the linear 
perturbation term we use 

AW n=1 = W e n Tf - W n=l = - + - , (39) 

w r 

and upon substituting (39) into (8) one clearly sees that full energy corrections are 
Ae n= i = Aen + Ae 12 = B — \ 2 /w 2 , together with the corresponding moderator 
function involving all the modifications 

A- = exp( ~ A ' >> ■ (40) 
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Hence, the first excited state wave function for (33) reads 

ipn=i = Xn=i0«=i = N exp (-wr 2 1 '4) exp (-Xr/w) , (41) 

from which it is clear that the present formalism produces exact results. However, 
the correction terms in KS approach [9] obtained for n — 1 state up to the second 
order given in (34) deviate slightly from the exact result (see Eq. (25) of [9]) as 
in the previous example, while it leads to the correct answer for the correction 
terms to the energy. 



3.4 infinite square well with linear perturbation (n=2 state) 

Finally, Kim and Sukhatme used their technique to investigate the n = 2 state 
of an infinite square well potential 

nr) - { I • m 

subject to a perturbation AV(r) = Ar + B. To calculate the energy level shifts 
and corrections to the unperturbed wave function 

Xn=2 = ^cos3r, (43) 
one first needs to define the corresponding superpotential 

W n=2 = = 3 tan 3r . (44) 

Xn=2 

To see once more the close relation between KS approach and the present for- 
malism, consider Eq. (11) in which AW 2 i can be obtained through (30) 

ATTr Asec 2 3r/vr 2 ,\ A / 1\ 
AW 21 = I -j-r 2 ) -- ^rtan3r + -J , (45) 

since the second order correction to the wave function in [9] appears in the 
form 



1 I 2 _ 1 



r tan 3r H — r 

4 / 9 



(46) 



Eq. (45) can also be checked out using (31) by setting the limits of the integral 
from — 7r/2 to 7r/2. Substitution of Eqs. (44) and (45) into (11) verifies the result 
Ae2i = B in [9]. Similarly, the second order energy level correction found in [9] 
is 

622 ~ 36 V 12 36 j ' ^ ' 
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which can in our theoretical framework be obtained by 

tt/2 

Ae 22 = - f - cos 2 3r AW^ (r) dr , (48) 

J IX 

-tt/2 

that is derived from (12). 

Again, the study of this example justifies that the formalism introduced in 
this letter for perturbative problems has a more general form, since the whole KS 
approach for the calculation of corrections due to perturbative terms is hidden 
in the subset of our formalism, namely the expansion of Eq. (8) leading to Eqs. 
(11) and (12) for the first and second order modifications reproduces same results 
as in KS formalism. Additionally, the mathematical treatment is much simpler 
in the frame of the present approach. 

4 Concluding Remarks 

No single approximation method available in the literature is ideal for every prob- 
lem. SSPT and LPT based theories avoids the Rayleigh-Schrodinger summation, 
but it can lead to nasty integrals and more effort in particular for excited states. 
The method is valuable when the integrals can be done exactly or by a reliable nu- 
merical procedure. Otherwise, the Rayleigh-Schrodinger summation, even when 
it does not give an exact answer, starts not to look so bad after all. This was 
the motivation behind the work introduced in this article. The present pertur- 
bation model appears in this respect to be superior for the excited states and 
provides a quick route to the calculation of all corrections within the frame of the 
perturbation theory, which considerably simplify one's calculational workload. 

Further, an attempt is made through this article to generalize SSPT and 
shown that the new formalism unifies SSPT- and LPT-based perturbation the- 
ories, which resolves difficulties and inefficiecies calculating in particular excited 
state corrections without having used tedious procedures. The power and el- 
egancy of the unified model which is, in a sense, complete are illustrated via 
specifally chosen four examples. We now have clear and explicit ways to get 
corrections to all energy levels and state wave functions for a given perturbed 
potential. 

As a concluding remark, SSQM has so far illuminated so many different areas 
of quantum mechanics that one might wonder what light perturbation theories 
based on supersymmetry might shed on the conventional techniques. In that 
hope, we are not about to be disappointed. We believe that this simple and in- 
tuitive method discussed through this paper would find a widespread application 
in the related area. The application of the method to other potentials involving 
Yukawa and Woods-Saxon potentials is in progress. 
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